function [dx, dy, dz] = mclm(x, mu, k)
    % MCLM模型:双Logistic映射 + 记忆阻器耦合
    
    diff = x(1) - x(2);
    dx = mu * x(1) * (1 - x(1)) + k * sin(x(3)) * diff;
    dy = mu * x(2) * (1 - x(2)) - k * sin(x(3)) * diff;
    dz = x(3) + diff;

end
